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ABSTRACT 





Context. The mysterious solar system object 133P/(7968) Elst-Pizarro is dynamically asteroidal, 
yet displays recurrent comet-Kke dust emission. Two scenarios were hypothesized to explain this 
unusual behavior: (1) 133P is a classical comet from the outer solar system that has evolved onto 
a main-belt orbit, or (2) 133P is a dynamically ordinary main-belt asteroid on which subsurface 
ice has recently been exposed. If (1) is correct, the expected rarity of a dynamical transition onto 
an asteroidal orbit implies that 133P could be alone in the main belt. In contrast, if (2) is correct, 
other icy main-belt objects should exist and could also exhibit cometary activity. 
Aims. Believing 133P to be a dynamically ordinary, yet icy main-belt asteroid, I set out to test 
the primary prediction of the hypothesis: that 133P-lik;e objects should be common and could be 
found by an appropriately designed observational survey. 

Methods. I conducted just such a survey — the Hawaii Trails Project — of selected main-belt 
asteroids in a search for objects displaying cometary activity. Optical observations were made of 
targets selected from among the Themis, Koronis, and Veritas asteroid families, the Karin asteroid 
cluster, and low-inclination, kilometer-scale outer-belt asteroids, using the Lulin 1.0 m. Small and 
Moderate Aperture Research Telescope System (SMARTS) 1.0 m. University of Hawaii 2.2 m. 
Southern Astrophysical Research (SOAR) 4.1 m, Gemini North 8.1 m, Subaru 8.2 m, and Keck I 
10 m telescopes. 

Results. 1 made 657 observations of 599 asteroids, discovering one active object now known 
as 176P/L1NEAR, leading to the identification of the new cometary class of main-belt comets. 
These results suggest that there could be ~100 currently active main-belt comets among low- 
inclination, kilometer-scale outer belt asteroids. Physically and statistically, main-belt comet ac- 
tivity is consistent with initiation by meter-sized impactors. The estimated rate of impacts and 
sizes of resulting active sites, however, imply that 133P-sized bodies should become significantly 
devolatilized over Gyr timescales, suggesting that 133P, and possibly the other MBCs as well, 
could be secondary, or even multigenerational, fragments from recent breakup events. 



Key words, comets: general - comets: individual: 1 33P/Elst-Pizarro - comets: individual: 
176P/L1NEAR - minor planets, asteroids - solar system: general 
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1. Introduction 

1.1. Motivation 

Comet 133P/(7968) Elst-Pizarro (hereafter 133P), an apparent Themis asteroid family member, 
has an asteroid-like Tisserand invariant of Ty = 3.16 and orbits in the outer main asteroid belt. 
Its cometary nature was first revealed on 1996 August 7 when a linear dust feature was observed 
trailing the object ( |Elst et ar||1996| l. This dust trail was found to be narrow, structureless, and over 
3 arcmin in length. No coma was observed around the nucleus. Previous images from 1979 and 
1985 ( |Marsden| [1996) [McNaught et a/.| [T996| and follow-up observations by |Offutt ef oT] ( [T997] l 
revealed a completely point-source-like nucleus with no apparent dust trail. 

The uniqueness of 133P's dust emission — no other asteroid had ever been seen exhibiting such 
behavior, and even 133P only demonstrated such behavior once — made it difficult to ascertain 
its true nature. Some believed the emission to be impact ejecta from a collision between 133P 



and another asteroid or group of asteroids (e.g., Toth 2000 1, while others believed it to be the 



result of the sublimation of volatile ices on 133P's surface (e.g., Boehnhardt et al. 1998 1. Impact 



scenarios were unable to plausibly explain the months-long duration of the dust emission episode 



( Boehnhardt et al. 1996 1, however, while sublimation required the apparently implausible scenario 
of surface ice surviving to the present-day on a main-belt asteroid in quantities sufficient for driving 
cometary activity. 

New insight was gained into 133P's nature when observations in 2002 revealed the return of 



long-lived dust emission ( Hsieh et al. 2004 1. Given the implausibility of impacts causing prolonged 
dust emission episodes on the same asteroid in the span of six years, these observations effectively 
ruled out hypotheses in which impacts were the sole cause of the observed dust emission. After 
considering other possible emission mechanisms and ruling out each in turn, Hsieh et al. concluded 
that 133P's dust ejection was most likely driven by the sublimation of volatile material, presumed 
to be water ice. 

This conclusion implied either that (1) 133P was a "lost comet", i.e., an object originally from 
the outer solar system that had evolved onto its current orbit via planetary encounters, dynamical 
resonances, or the non-gravitational influence of asymmetric cometary mass loss, or (2) 133P was 
an "icy asteroid", i.e., a dynamically ordinary, native member of the asteroid belt on which pre- 
served, buried ice had been recently exposed, perhaps excavated by an impact. If the "lost comet" 
hypothesis was correct, the low likelihood of a comet undergoing a dynamical transition from the 



outer solar system onto a 133P-like orbit (e.g., Ipatov & Hahn 1997 Fernandez et al. 2002| l could 



* Some of the data presented herein were obtained at the W. M. Keck Observatory, the Gemini Observatory, 
Subaru Telescope, National Optical Astronomy Observatory (NOAO) facilities at the Cerro Tololo Inter- 
American Observatory, and Lulin Observatory. Keck is operated as a scientific partnership among the 
California Institute of Technology, the University of California, and the National Aeronautics and Space 
Administration, and was made possible by the generous financial support of the W. M. Keck Foundation. 
Gemini is operated by the Association of Universities for Research in Astronomy, Inc., under a coopera- 
tive agreement with the National Science Foundation (NSF) on behalf of the Gemini partnership. Subaru 
is operated by the National Astronomical Observatory of lapan. NOAO and Cerro Tololo are operated by 
the Association of Universities for Research in Astronomy, Inc., under co-operative agreement with the 
NSF. Lulin is supported and was made possible by the National Science Council of Taiwan, the Ministry 
of Education of Taiwan, and National Central University. 
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mean that such objects are rare or even non-existent, and that 133P could be the only object of its 
kind, as was observed. If the "icy asteroid" hypothesis was correct, 133P would not be expected 
to be unique. Other main-belt asteroids could also possess subsurface ice that could begin subli- 
mating if exposed by impacts. In this case, 133P's uniqueness could simply be due to the difficulty 
of discovering such objects serendipitously. An appropriately targeted search could counter this 



difficulty ( Hsieh & Jewitt 2006a I. Believing that this second hypothesis could be correct, I set out 
to conduct just such a targeted search. In this paper, I describe the design and results of that survey, 
the Hawaii Trails Project (HTP). 



1.2. Background 



Various evidence has been found for past and even present water in main-belt asteroids {cf. Hsieh 
|& JewittI |2006a[ and references within), but this fact has not meant that active dust emission due 
to the sublimation of water ice is commonly observed. Hsieh & Jewitt discuss a number of pos- 
sible reasons for this inconsistency. Perhaps simply too few 133P-like objects had been observed 
sufficiently deeply and sufficiently often. CoUisional excavation of subsurface ice reservoirs could 
occur only rarely. Even when collisions do occur, sublimation may not necessarily follow due to 
the probable non-uniform distribution of subsurface ice reservoirs, and even the successful trig- 
gering of sublimation could result in dust emission that is simply too weak to detect from Earth. 
Finally, observable dust emission is likely to be intermittent and also have a finite lifetime, limited 



by either local devolatilization of the active site or mantling (Jewitt 1996 2002 and references 
within), farther reducing the probability of its detection. 

In addition to the above considerations, there are indications that detectable activated asteroids 
may occupy a narrow range of sizes. The rate of impact excavations of an asteroid depends on 
its collisional cross-section: large bodies will be struck, and thus activated, more often than small 
bodies. Volatile material should also survive longer against solar heating on larger asteroids because 
it can be buried at greater depths. However, Hsieh et al. ( 2004| l found that the low dust ejection 



velocity (vgj ~ 1 - 2 m s" ) for 133P was extremely similar to the gravitational escape velocity 
(Vesc ~ 1.5 m s"'), suggesting comparable dust emission would only be observable on bodies of 
similar or smaller sizes. Comparable emission on larger bodies would be unable to escape the 
gravity of those bodies and thus never become observable. This scenario may explain why Ceres 
is not observed to emit dust, despite indications of surface ice ( [Lebofsky et al.\ |1981[ |Vernazza~ef| 
al. 2005 1 and even possible water vapor emission ( A'Hearn & Feldnian]|1992 1. Despite providing 



less insulation from the Sun than larger bodies, smaller bodies are also less effectively heated from 
within by ^''Al (since they radiate that heat away more efficiently). In summary, to exhibit dust 
emission, it seems that icy asteroids need to be large enough to preserve subsurface ice against 
solar heating and be struck relatively frequently by impactors, but not so large that weak dust 
emission cannot escape gravity or ice does not survive early radioactive heating. Such conflicts in 
size preferences suggests that 133P-like asteroids may be even rarer and more difficult to discover 
than initial considerations might indicate. 
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1.3. Survey Design 

The distribution of present-day ice in the asteroid belt is not well understood. The so-called snow 
line refers to the distance from the Sun at which the temperature of the protosolar disk was below 
the condensation temperature of water, allowing ice to become incorporated into accreting plan- 
etesimals. Observations of asteroids suggest that the snow line likely existed ~2.5 AU from the 



Sun ( Gradie & Tedesco 1982 Jones ef a/. 1990 1, though theoretical models have suggested it may 



have been as close as the orbit of Mars ( |Sasselov & Lecar[ [2000 [ [Ciesla & Cuzzi| |2006[ |Lecar et\ 
al. 2006| l, suggesting objects throughout the asteroid belt could have formed icy. Much of this ice 
has since undergone thermal processing, producing the hydrated minerals now observed on many 
main-belt asteroids (see |Rivkin et al.\\2002\ . Detections of such minerals decline in the outer belt 
beyond ~3 AU (notably, where 133P is found), suggesting that pristine, unaltered ice could still 



exist in distant main-belt objects (Jones ef a/. 1990 Scott &Krot 2005 1. This hypothesis suggests 
that a search for 133P-like icy asteroids would be best served by focusing on objects with semi- 
major axes of a > 3 AU, as the HTP does indeed do, although objects with a < 3 AU are also 
considered. 

Given the environmental and physical conditions discussed above that are suspected of giving 
rise to cometary activity in 133P, in the search for objects exhibiting similar activity, the following 
target categories were considered: 

- Themis family asteroids: The Themis family is likely the result of the breakup of a parent 



asteroid about 400 km in diameter ~2.5 Gyr ago (Marzari et al. 1995 Tanga et al. 1999 



Nesvorny et al. 2003[ l and is one of the largest and most statistically robust asteroid fami- 



lies known {e.g., Carusi & Valsecchi 1982[ Zappala et al. 19901. Given their origin from a 
common parent, Themis family members are thought to be compositionally homogeneous, as 
corroborated by observational studies showing that the family is dominated by primitive C-type 



asteroids that also exhibit signs of aqueous alteration ( Bell 1989 Florczak ef aZ. 1999 Ivezic 



et al. 2002 Mothe-Diniz et al. 2005 [ ). The family is also characterized by a high rate of col- 



lisions relative to the general main-belt population (Farinella & Davis 1992[|DeirOro et al. 



|2001| l. Given 133P's dynamical and possible compositional similarity to Themis family mem- 
bers and the high frequency of collisions within the family, the Themis family was judged to 
be a natural place to look for 133P-like asteroids. Target objects were selected by conducting a 
hierarchical clustering analysis (Zappala et al. 1990 1994| l on proper elements obtained from 
the AstDys website (http://hamilton.dm.unipi.it/astdys/ ; Knezevic & Milani 2003| l and using 
a cutoff value of 5v' = 70 m s~' with respect to (24) Themis. 
- Low-inclination outer-belt asteroids: The identification of the main-belt comet (MBC) class 
(Hsieh & Jewitt 2006b| l during the course of the HTP led to a re-evaluation of the survey's 
target selection criteria. In particular, the semimajor axis and inclination of MBC P/2005 
Ul (Read) (hereafter P/Read) closely matched those of the Themis family, but its eccentric- 
ity did not. This discovery led to the hypothesis that 133P's low-inclination, and not neces- 
sarily its specific association with the Themis family, could be primarily responsible for its 
activity. Collision rates are in fact known to be enhanced at low inclinations ( [Farinella &| 
Davis 1992[ l. Targets were selected by applying limits (3.0 AU < a < 3.3 AU, e < 0.3, 



and ; < 3.0°) on osculating orbital elements obtained online from Lowell Observatory 
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( |ftp://ftp.lowell.edu/pub/elgb/astorb.dat| , permitting the inclusion of even recently discovered 
asteroids for which proper elements were not yet known. Limits were chosen to span the orbital 
element distribution of the first three known MBCs, and in fact, also completely encompasses 
the Themis family. 

- Koronis family asteroids: As another populous main-belt asteroid family characterized by 



high collisions rates (Farinella & Davis 19921, the Koronis family was judged to be another 
plausible place to look for 133P analogs. Like the Themis family, the Koronis family is the 



result of the catastrophic disruption of a large parent body ~2 Gyr ago (Binzel 1988 Marzari 



et al. 1995| l. Unlike the Themis family, however, the Koronis family is dominated by S-type 



asteroids ( |Bell| |1989[ [Binzel et oT] |1993[ [Mothe-Diniz et al.\ |2005[ ), which along with their 
meteoritic analogs, the ordinary chondrites, do not exhibit spectroscopic evidence of signifi- 



cant aqueous alteration {e.g., Rivkin et al. 20021. Nonetheless, due to the high collision rates 
within this family and the possibility that these asteroids may not actually be completely an- 
hydrous (Grossman et al. 2000[ Keil 20001, Koronis family asteroids were still judged to be 
viable 133P-analog candidates. Koronis family targets were selected by a hierarchical cluster- 
ing analysis, using a cutoff value of 6v' = 50 m s ' with respect to (158) Koronis. 
- Karin cluster asteroids: The Karin cluster is contained within the Koronis family and is 



likely the result of a recent fragmentation event ~5.8 Myr ago (Nesvorny et al. 2002 20061. 
Spectroscopic studies have since affirmed both the compositional homogeneity and the young 



age of the cluster ( Sasaki ef a/. 2004 Brunetto a/. 2006 Vernazza e/ a/. 2006 1. In the con- 
text of this survey, the Karin cluster was interesting because if its parent body only recently 
broke apart, ice that had been preserved deep inside could suddenly be located at a much shal- 
lower depths on the resulting fragments. Even relatively weak subsequent impacts could then 
easily initiate sublimation. Like the rest of the Koronis family, however, the Karin cluster is 
dominated by S-type asteroids and so is not an ideal sample from a mineralogical point of 
view for searching for sublimation-driven dust emission. Given the likelihood of particularly 
young surfaces, though, these objects were judged to be worthy of at least cursory attention. 
Karin cluster targets were selected by a hierarchical clustering analysis, using a cutoff value of 
dV = 10 m s ' with respect to (832) Karin. 



- Veritas family asteroids: Like the Karin cluster, the Veritas family, first identified by Zappala 



et al. ( 1995| l, is likely the result of a recent breakup in the asteroid belt, perhaps just ~8.5 Myr 



ago ( Nesvorny et al. 2003 [ Tsiganis et al. 2007 1. Intriguingly, such a breakup corresponds with 
a spike in ^He concentrations found in sea-floor sediment records from that time period ( Farley 



et al. 2006| l. The young age of the Veritas family made it interesting for this survey for the 



same reasons as the Karin cluster. Unlike the Karin cluster, however, the Veritas family appears 
to be dominated by C-type asteroids ( |Mothe-Diniz et oL 2005 i, the same type of asteroids 
that dominate the Themis family, making it a particularly promising place to search for 133P 
analogs. Veritas family targets were selected by a hierarchical clustering analysis, using a cutoff 
value of 6v' = 40 m s"' with respect to (490) Veritas. 



Targets from these categories were selected for each observing run based on visibility, expected 
visual magnitude (18 ;$ $ 20 mag targets for UH 2.2 m runs, brighter targets for smaller 
telescopes or in non-photometric conditions, and fainter targets for larger telescopes), distance from 
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the moon, and distance from the galactic plane. I also sought to observe asteroids matching 133P 
closely in physical size (as parameterized by absolute magnitude, Hy), or approximately kilometer- 



scale, given the issues discussed m{\.2 The aim was to achieve moderately deep imaging for as 
many asteroids as possible. As such, exposure times were selected such that signal-to-noise ratios 
of at least S /N = 20 for 93.5% (614 out of 657) of the surveyed targets, and at least S /N ^ 10 for 
99.2% (652 out of 657) of the surveyed targets were reached. 

While most observations were made of unique asteroids, repeat observations of individual as- 
teroids were not expressly avoided. Given the expected transient nature of any cometary activity, 
continued monitoring of comet candidates is in fact, in principle, vital to fully characterizing their 
active natures. In a larger survey, all objects would have been observed multiple times, but for the 
relatively small-scale HTP, repeat observations were made only occasionally. 



2. Observations 

Several ground-based optical observatories were used to conduct this survey: the University of 
Hawaii (UH) 2.2 meter telescope, the 8.1 m Gemini North Observatory, the 8.2 m Subaru tele- 
scope, and the 10 m Keck I Observatory on Mauna Kea in Hawaii, the 1.0 m telescope operated by 
the Small and Moderate Aperture Research Telescope System (SMARTS) consortium at the Cerro 
Tololo Inter- American Observatory (CTIO) and Cerro Pachon's 4.1 m Southern Astrophysical 
Research (SOAR) telescope in Chile, and the Lulin 1.0 m telescope in Taiwan (Table [TJ. Data 
was obtained on 88 separate nights (Table |2]l. 

Either a Tektronix 2048x2048 pixel charge-coupled device (CCD) or the Orthogonal Parallel 
Transfer Imaging Camera (OPTIC) was used at the f/10 focus of the UH 2.2 m telescope. Other 
instruments used on Mauna Kea telescopes were the Gemini Multi-Object Spectrograph (GMOS) 
( Hook et al. 2004| l in imaging mode on Gemini, the Subaru Prime Focus Camera (Suprime-Cam) 



(Miyazaki et al. 2002[ l on Subaru, and the Low-Resolution Imaging Spectrometer (LRIS) (Oke 



et al. 1995 1 on Keck. Observations using the CTIO 1.0 m telescope employed either an Apogee 



512x512 CCD or an STA 4064x4064 CCD (Y4KCam), while observations on the SOAR telescope 



employed the SOAR Optical Imager (SOI) ( Schwarz et al. 2004 1. Observations using the Lulin 



1.0 m telescope were made using a Vers Array :1300B CCD (Kinoshita et al. 2005 i. All observa- 
tions were obtained through broadband filters approximating the Kron-Cousins BVRI photometric 
system except for observations at Gemini which employed g'r'i'z' broadband filters similar to those 
used for the Sloan Digital Sky Survey (SDSS). Further telescope details are listed in Table[T] 

Standard image preparation (bias subtraction and flat-field reduction) was performed on all 
data. Flat fields were constructed either from dithered images taken nightly of the twilight sky or 
images of the illuminated interior of the telescope dome. Photometry on standard stars and target 
objects was performed by measuring net fluxes contained within circular apertures with optimum 
sizes determined from curve of growth analysis and dependent on the nightly seeing. The sky 
background to be subtracted was determined from the median pixel value within a circular annulus 



surrounding each central aperture. Photometry was calibrated to Landolt ( 1992 1 standard stars to 
produce absolute photometry for the target objects. Photometric uncertainties are estimated to be 
0.1 mag for objects observed in photometric conditions and 0.5 mag for objects observed in non- 
photometric conditions (see Table|2]i. 



H. H. Hsieh: The Hawaii Trails Project 7 

3. Analysis and Results 

3.1. Activity Assessment 

To search for low-level cometary activity, all images for a given target on a single night were aligned 
on the object's photocenter using linear interpolation and then summed to produce a single compos- 
ite image. Asteroids with possible cometary activity were identified visually from these composite 
images, and then the corresponding individual images were inspected to assess the plausibility of 
the detection (e.g., whether the morphology of a suspected cometary feature matched that expected 
for comets, whether the feature consistently appeared in each individual image, whether the motion 
of the feature and object against the background star field were identical, and whether there were 
any signs that the features could be attributed to scattered light or imperfect flatfielding). After this 
close inspection, if an object showed credible evidence of being active, confirmation of activity 
was then attempted using deep follow-up imaging. Observations for which activity was impossi- 
ble to assess and were therefore unsuitable for analysis (e.g. those where the object was strongly 
affected by scattered light from a nearby bright star, or where poor telescope stability resulted in 
inconsistent point-spread functions among individual images of the same object) were discarded. 
Various systematic, quantitative means of searching for cometary activity in HTP data were 



considered, including one-dimensional surface brightness profile analysis (e.g., Luu & Jewitt 1992 



Hsieh & Jewitt 20051, radial surface brightness profile analysis (e.g., Luu & Jewitt 1990 1, and 



trail searching (Hsieh & Jewitt 2005 1 ). One-dimensional surface brightness profile analysis can be 
useful in that it accounts for trailing of point source templates (i.e., field stars) caused by non- 
sidereal tracking and long exposure times of each target image. Trailing effects are avoided by 
measuring the surface brightness profiles of targets and field stars in the direction perpendicular 
to the apparent motion of the target. This technique is best-suited for detecting coma that extends 
radially in all directions from an object, or directed emission not aligned with the direction of 
the object's apparent motion, but would miss directed emission oriented along the direction of an 
object's apparent motion (as was observed for 133P). 

Radial surface brightness profile analysis is sensitive to cometary features extended in any di- 
rection, but requires the use of circularly-symmetric point source templates. Field stars from target 
images were generally unsuitable due to trailing, requiring the use of either standard star images 
or separately-observed sidereally-tracked target images to obtain reference star templates, but both 
of these alternatives were problematic as well. Even with the use of tip-tilt guiding, continuous at- 
mospheric fluctuations and small but non-negligible mechanical instabilities in the telescopes used 
(particularly while tracking non-sidereally) resulted in systematically wider point-spread functions 
(PSFs) for long-exposure target images than short-exposure standard star images. Those same im- 
age quality fluctuations also adversely affected the utility of separate sidereally-tracked target im- 
ages by making it difficult to distinguish ordinary image-to-image PSF fluctuations from the small 
deviations that could actually indicate the presence of activity. This technique also averages any 
dust emission over all radial directions, reducing sensitivity to directed emission. 

Searching for dust emission specifically like 133P's dust trail, which could form if an object 



emits large (^ 10 fim) dust particles at low velocities (~ 1 - 2 m s ) as 133P did (Hsieh et al. 



2004[ l, can be done by searching for excess flux above sky background along the object's orbital 



plane (the expected direction of a trail). This technique is not sensitive to near-nucleus activity. 
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however, and furthermore, is highly susceptible to contamination from field stars or galaxies due 
to the large amount of sky that must be examined for evidence of trails. 

All three techniques are prone to false detections, which can arise from nearby faint field stars 
or galaxies masquerading as active features, contamination of sky background samples, imperfect 
flatfielding, scattered light from nearby bright or saturated field stars, and PSF variability between 
images or even within the same image (i.e., as a function of position on the CCD). Thus, any auto- 
mated detection scheme still requires confirmation by visual inspection. Given the focused nature 
of this survey (targeted deep imaging of single objects), the resulting small sample size (~650 
composite images to be evaluated over 3 years), and the significant human intervention needed 
to supplement any of the quantitative detection algorithms considered, I opted to simply visually 
inspect each target for evidence of activity. This approach is certainly less feasible for searching 
for comets in much larger data sets, particularly the ones that will be produced by upcoming sur- 
vey telescopes like the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; 
[Hodapp et a/.[|2004| and the Large Synoptic Survey Telescope (LSST; |Jones et a/.||2009| l. For the 
HTP, however, direct visual inspection was judged to be the most efficient (given the small sample 
size) and effective method for evaluating the active nature of the targets in this survey. 



3.2. Survey Results 

In total, 599 unique asteroids were observed as part of the HTP, with 544 asteroids observed once, 
52 asteroids observed twice, and 3 asteroids observed three times, totaling 657 separate observa- 
tions (Table[3]). Hereafter, all observations (including those of objects observed multiple times) are 
considered independent of one another. This approach is taken since it is not known where along 
their orbits other 133P-like asteroids could exhibit activity, and also since inactivity observed for 
an asteroid at one time does not necessarily rule out activity at another time. 

The numbers of targets observed from each target category are shown in Table |4] and Fig.[T[ 
for which it should be noted that all Themis family members are also counted as low-inclination 
outer belt asteroids, and all Karin cluster members are also counted as Koronis family members. 
The Veritas family is independent of all other target categories. Slightly more than half (389 of 
657) of observed targets were Themis family members, reflecting the prediction that these were the 
objects among which other 133P-like asteroids would most likely be found. This histogram, and 
others following it, also show the relative numbers of targets for which shallow, medium, and deep 
surface brightness detection limits were reached. 

The semimajor axis versus eccentricity distribution and semimajor axis versus inclination dis- 
tribution of surveyed targets in the context of the background main-belt population are shown in 
Figs. |2]and[3] respectively. These distributions primarily reflect the family association or orbital 
element criteria of HTP target categories ( ^1.3[ l. 



3.3. Discovery of the Main-Beit Comets 



The most significant result from the HTP was the discovery of cometary activity in asteroid 11 8401 
(1999 RE70; subsequently designated 176P/LINEAR; hereafter 176P), on 2005 November 26 by 
Gemini North (Hsieh & Jewitt 2006b Hsieh et al. 2006[ l. Coincidentally, this discovery was 
made just a month following the serendipitous discovery of P/Read ( Read et al. 2005 1, another 



H. H. Hsieh: The Hawaii Trails Project 9 

cometary object orbiting completely within the main asteroid belt, on 2005 October 24 by the 
0.9 m Spacewatch telescope. Together, 133P, P/Read, and 176? formed a new class of objects 
which were dubbed MBCs ( |Hsieh & Jewittj |2006b| ). Following the completion of the HTP, a 
fourth MBC, P/2008 Rl (Garradd) (hereafter, P/Garradd; [Garradd et aT\ [2008] ), was discovered 
on 2008 September 2 by the 0.5-m Uppsala Schmidt telescope at Siding Spring. The properties of 
the currently-known MBCs are listed in Table |5] All other HTP observations (including a previ- 
ous observation of 176P using the Lulin 1.0 m telescope one month prior to the discovery of the 
object's cometary nature by Gemini) were judged to show no evidence of cometary activity. 

The extremely limited data set used to discover the currently-known population (three members 
were discovered serendipitously by 0.5 m to 1.0 m telescopes, and another was discovered by a tar- 
geted survey composed of only 657 asteroid observations) impUes that a much larger undiscovered 
population exists. The existence of multiple known MBCs and the implication of a much larger 
as-yet undiscovered population runs counter to the prediction made by the "lost comet" hypothesis 
that 133P should be alone in the main belt, making the "icy asteroid" hypothesis, which predicts 
that such objects should be widespread in the main belt, much more likely (§ 1. 1 1. 



4. Discussion 

4.1. Selection Effects and Biases 

It should be emphasized that the HTP was not designed as a random sampling of the main asteroid 
belt and cannot be interpreted as such. On the contrary, it was designed to maximize the probability 
of discovering 133P-like asteroids and so employed very specific physical and dynamical target 
selection criteria. The HTP was also subject to a number of ordinary observational biases. 

Histograms showing the distributions of various physical and dynamical properties and obser- 
vational circumstances for surveyed targets are shown in Figs.|4]and[5] For reference, the locations 
of the known MBCs (for Fig. |4| or the range of locations of the known MBCs while they were 
observed to be active (for Fig. [5]l in these distributions are marked. These histograms reflect the 
following intentional selection effects or observational biases: 

- The distributions of semimajor axes (Fig.|4^), eccentricities (Fig.|4]3), and inclinations (Fig. [4];) 
of surveyed targets strongly reflect the focus on specific asteroid families, particularly the low- 
inclination, moderate eccentricity Themis family at ~ 3.1 - 3.2 AU and the low-inclination, 
low-eccentricity Koronis family at ~2.9 AU (also see Figs. |2] and [3]l. 

- The distribution of absolute magnitudes (Fig.HH) of surveyed targets reflects the focus on aster- 



oids similar to 133P (Hy = 15.9 mag; Hsieh et al. 2009a I. It is also affected, however, by the 



fact that smaller targets are fainter and are more difficult to image to the same depth as larger, 
brighter targets. The known population of small objects is also quite incomplete, again due 
to the faintness of these objects and the difficulty of discovering them. Larger asteroids were 
furthermore more likely to have well-known orbits and calculated proper orbital elements, nec- 
essary for computing family linkages on which initial target selection relied. 
- The increased difficulty of observing fainter objects farther from perihelion (v - 0°) meant that 
fewer were observed, as can be seen from the heliocentric distance and true anomaly distri- 
butions of HTP observations (Figs.|5^ andjSj)). The resulting oversampling of near-perihelion 



objects means that even if the seasonal activity modulation hypothesis discussed above (^-2 1 is 
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correct and there is an equal likelihood of activity along an asteroid's entire orbit, active objects 
are still more likely to be near perihelion when they are observed. For this reason, the dis- 
covery of near-perihelion activity in this survey should not be interpreted as confirmation that 
activity in 133P-like asteroids is correlated to perihelion passages. Similarly, the correlation of 
activity of serendipitously discovered MBCs to perihelion passages should not necessarily be 
interpreted as real either since this correlation may simply be due to the fact that the objects 
are closer to observers near perihelion and their activity is therefore brighter and more easily 
detected. 

- The low inclinations of most targets (Fig. [4|;) necessarily meant that they were also observed 
at low orbit plane angles (Fig. |5};). This trend is significant because low orbit plane angles 
are where dust emission is most easily observed. When dispersed by solar radiation pressure, 
emitted dust will tend to spread out in a plane behind an object and aligned with its orbit. Thus, 
when the Earth is also aligned with the object's orbit plane, the optical depth of the emitted dust 
is maximized. Indeed, when observed to be actively emitting dust, the first three known MBCs 
were at low orbit plane angles with respect to the Earth (la^/l < 0.6°). 

- The distribution of apparent magnitudes (Fig. |5}l) primarily reflects the fact that the majority 
of HTP observations were made using the UH 2.2 m telescope and objects with visual R-hand 
magnitudes of 18 ^ mu ^ 20 mag were judged to be the most effectively and efficiently 
observed with that telescope aperture. Fainter asteroids were observed with larger telescopes 
such as Subaru, Gemini, and Keck I, but these targets typically were not imaged as deeply as 
brighter objects due to the desire to also observe as many of these smaller asteroids as possible. 



It should be noted that P/Read was serendipitously discovered and so its strong orbital simi- 
larities to 133P and 176P cannot be easily explained by selection effects or observational biases. 
Its discovery may in fact indicate a true predisposition of the region near the Themis family for 
harboring MBCs. It should also be noted, however, that the discovery of P/Garradd in a completely 
unexpected, and unsurveyed, region of the main belt underscores the point that the HTP cannot and 
should not be viewed as a representative sampling of the cometary content of the entire main belt. 

Since the discovery of 176P's cometary nature as part of the targeted HTP, at least two untar- 
geted surveys have been conducted to search for other MBCs. [Gilbert & Wiegert| ( |2009| l searched 
12390 main-belt objects observed as part of the Very Wide segment of the Canada-France-Hawaii 
Telescope (CFHT) Legacy Survey, finding one known comet (active Centaur 166P/NEAT) and one 
unknown comet with an unknown orbit, but no active objects with confirmed main-belt orbits. For a 
smafl portion of their survey (952 targets), automated techniques, some similar to those considered 
for the HTP, were used to search for cometary activity. As in the HTP, though, the majority of their 
targets (11438 targets) were simply assessed visually. Both cometary objects were found by visual 
inspection. From these results, they estimated an upper limit of ~ 800 weakly active MBCs among 
objects larger than 1.5 km in diameter with semimajor axes between 3.0 AU and 5.0 AU, consistent 



with the results found by the HTP (54.3 i. Meanwhile, Sonnett et al. (2008i searched 1850 objects 



observed with the CFHT as part of the Thousand Asteroid Light Curve Survey (Masiero et al. 



2009| l. Several dififerent automated techniques were employed to search for activity, but ultimately 



no evidence of activity was found in any of those objects. 
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Such untargeted surveys are clearly better suited than the HTP for developing generalizations 
about the prevalence of MBCs throughout the asteroid belt in that they do not focus on specific re- 
gions of orbital parameter space. The correlated disadvantage of the untargeted approach, however, 
that the probability of discovering MBCs is much lower than that of targeted surveys like the HTP, 
as is reflected by the lack of detections in either untargeted survey, despite both having much larger 
target sample sizes than the HTP. 

4.2. Survey Sensitivity 

The image depth to which candidate objects should be observed was a key uncertainty at the 
start of this survey. When an active 133P was observed on 2002 November 5, the surface bright- 
ness of a representative section of the trail (just outside the seeing disk of the nucleus) was 
S = 25.2 mag arcsec"^. In previous observing runs that fall, 133P's trail was brighter Thus, its 
2002 November 5 surface brightness could be considered a minimum surface brightness sensitiv- 
ity for future observations searching for similar activity. This need for deep imaging had to be 
balanced with maximizing the target sample size to maximize the likelihood of finding an object 
that was actually active. In the end, ~75% (488) of the 657 observations made had 3cr surface 
brightness detection limits of 2;„„ = 25.2 mag arcsec"^ or better, though many observations went 
even deeper, with more than 25% (170 objects) having 3cr surface brightness detection limits of 
I,iim > 26.0 mag arcsec"^ (Table [3]l. 

The discovery of 176P's cometary nature during this survey was made with two 120 s exposures 
(which, stacked together, had a 3cr surface brightness detection limit of 2;,„, = 26.8 mag arcsec"^) 
on 2005 November 26 on the 8 m Gemini telescope. In these observations, the tail had a typical sur- 
face brightness ofL- 25.3 mag arcsec"^, a surface brightness which should have been detectable 
by ~70% (466) of HTP observations. One cautionary note is that 176P was also observed on 2005 
October 24 using the Lulin 1.0 m telescope just one month before the discovery of its cometary 
activity by Gemini on 2005 November 26. No activity was detected in the Lulin data, which had 
a surface brightness detection limit of S/,„, - 25.3 mag arcsec"^, indicating that 176P's tail could 
simply have been fainter then, or perhaps poor seeing (~ r.'5 for the Lulin data, compared to ~ 0'.'9 
for the Gemini data) could be to blame for the lack of observed activity. Only ~75% (484) of HTP 
observations were made on nights when the seeing was < 1'.' 1, and of these observations, 427 (65% 
of the total survey) had S/„„ > 25.3 mag arcsec"^. 

To provide a physical basis of comparison, the ratio of dust surface area to the total nucleus 
scattering cross-section, Cd/C„, contained within 1x10^ km at the geocentric distance of each target 
object, that is represented by 2/,m is also calculated and shown in Table[3] This parameter is roughly 
correlated with S /N and is computed to relate the activity detection limit of an observation (which 
is largely dependent on the total effective exposure time, size of telescope used, and proximity 
of the Moon or bright stars, but not the target itself) to the size and brightness of the target being 
observed. For reference, CdlC„ - 0.0039 for 133P's trail when observed by the UH 2.2 m telescope 
on 2002 November 5, while Cd/C„ = 0.0022 for 176P's tail when observed by Gemini on 2005 
November 26. Thus, activity as intense as 176P's activity, relative to the size and brightness of the 
nucleus, should have been detectable in ~70% of HTP observations (where 466 of 657 observations 
had surface brightness detection limits equivalent to Cd/C„ < 0.022), while activity as intense as 
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133P's activity, relative to its nucleus, should have been detectable in ~9Q% (588 of 657) of HTP 
observations. 



4.3. MBC Population Size 

A natural question to ask following the completion of the HTP is how many other active MBCs are 
impUed to exist? This question is rather tricky to answer given the various ways that HTP results 
can be scaled to the general main-belt population and the difficulty of interpreting the statistical 
implications of a single detection. To start out, however, the probability, Pam, that a randomly 
selected object in a given population will be active at a given time can be estimated by 

Pactv — ■ forb (1) 

"for 

where riactv is the number of active objects in the population, n,ot is the total number of objects in 
the population, and f„rb is the fraction of an object's orbit during which it is actually active. For this 



analysis, 1 adopt the active orbit fraction of 133P of forb ~ 0.25 (Hsieh et al. 2004 and references 
within) to be typical of all MBCs. 

Observationally, Pactv can be estimated from 

Pacn. = — (2) 

where ridet is the number of detections of activity in a particular sample set, and rtsun- is the total 
number of candidate objects observed that belong to that sample set. Thus, combining Equations[T] 
and|2] one obtains 

^survjorb 

The HTP detected one active MBC {ndet - 1), with the survey size ranging from Hsurv - 657, 
if all surveyed objects are included, to Usun- - 427, if only observations judged to be sufficiently 
sensitive and made in sufficiently good seeing conditions to have permitted the detection of 176P- 
like activity ( ^4.2| i are considered. Hereafter, only the latter group of objects, where risurv - 427, is 
considered. 

As of 2009 April 1, there are approximately 450,000 known main-belt asteroids. If the entire 
asteroid belt is considered to be the candidate population and using forb ~ 0.25 (assuming 133P to 
be typical of all MBCs), Equation [3] indicates that there could be ~4200 MBCs in the main belt. 
This result, however, suggests that all main-belt asteroids are equally likely to be MBCs, which is 
likely not the case. 

Refining this analysis, one might note that, just like 133P, 176P was discovered among the 
highly collisionally-active, C-type-dominated Themis family, and therefore Themis family targets 
may be a more appropriate sample from which to scale these results. In this case, «„„ = 2271 
and risurv - 262, giving nam ~ 35. Like scaling to the entire main-belt population, this calcu- 
lation assumes that all Themis family members are equally likely to display activity. As before, 
however, this is likely not the case ( ^1.2| i: small asteroids are more likely to display observable 
evidence of sublimation due to their small escape velocities and higher likelihood of being recently 
collisionally-produced fragments ( |Cheng| [2004| l that may have near-surface ice deposits. Thus far, 
this hypothesis is supported by the evidence. All four currently-known MBCs have Hy > 15. If 
only Themis family asteroids in this size range are considered, using nsurviHy > 15.0) = 205 and 
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fitoAHv ^ 15.0) = 358, one finds that there may only be ~10 MBCs that are detectable via the type 
of optical observations that were done in this survey. 

Refining this analysis yet again, one might note that unlike 133P and 176P, P/Read has an oscu- 
lating eccentricity that places it beyond the approximate bounds of the Themis family. It is possible 
that when proper elements become available for P/Read, it will in fact be found to be associated 
with the Themis family. Nonetheless, immediately following P/Read's discovery, 1 began survey- 
ing other asteroids that, hke P/Read, matched the Themis family in semimajor axis and inclination, 
but not in eccentricity. These are the HTP's low-inclination outer-belt targets. Using this sample 
(which only includes objects with Hy > 14.0) where n,ot - 8647 and nsurv = 327 (the total number 
of observed low-inclination targets with Hy > 14.0), I estimate that nactv ~ 100 among kilometer- 
scale, low-inclination, outer-belt asteroids. Using analogous lines of reasoning for the other target 
categories surveyed as part of the HTP where rtdet < 1 and only targets with Hy > 14.0 are con- 
sidered, 1 estimate that riacn < 150 for the Koronis family, riacn- < 40 for the Karin cluster, and 
riacrv < 15 for the Veritas family. 

The design of the HTP as a highly biased, targeted survey means that it cannot be used to 
infer the prevalence of MBCs in the general main belt population. There are large regions of the 
main belt that were not surveyed (c/.. Figs. |2] and [3]) and the recent discovery of the fourth-known 
JVIBC, P/Garradd, shows that MBCs can certainly exist in those unsurveyed regions (see Fig. |4]). 
Furthermore, preliminary results suggest that the nuclei of P/Read and P/Garradd are both ex- 



tremely small (sub-kilometer-scale; Hsieh et al. 2009b Jewitt et al. 2009 1, suggesting that, if not 
for their prodigious cometary activity, they would never have been discovered due to their extreme 
faintness. These results further suggest that more as-yet unknown sub-kilometer-scale cometary as- 
teroids might be discovered in the future, but without better knowledge of the size of the main-belt 
population at this size scale and a systematic, statistically significant survey of this population, it is 
difficult to predict how many of these small MBCs there might be. 



Jewitt et al. (2009 1 found that while P/Garradd is characterized by an asteroidal = 3.217, its 



orbit is unstable on a timescale of ~20 Myr, due largely to the influence of the nearby 8:3 mean- 
motion resonance (2.706 AU from the Sun) with Jupiter and the ve secular resonance (near 2.1 AU 
from the Sun). Because of these two facts, they hypothesize that the object may have migrated to its 
current location from a source region in the outer asteroid belt, perhaps where the other three known 



MBCs currently reside. This hypothesis is slightly problematic as de EHa & Brunini ( 2007 1 have 
actually found that mixing should be minimal among zones of the main asteroid belt dehneated by 
the strongly chaotic 3:1 (2.5 AU from the Sun), 5:2 (2.8 AU from the Sun), and 2:1 (3.27 AU from 
the Sun) mean motion resonances with Jupiter. In a study of V-type asteroids in the middle main 



belt, however, Roig et al. (2008 1 found that a non-negligible number of asteroids could be capable 
of migrating outward across Jupiter's 3:1 resonance, and thus it is conceivable that P/Garradd could 
be an icy outer-belt asteroid that analogously migrated inward across Jupiter's 5:2 resonance. 

From these results, it is clear that the small number of currently-known MBCs greatly hinders 
any attempt to estimate just how common they are in the main asteroid belt, and that better con- 
straints may not become available until deep-imaging, all-sky survey telescopes like Pan-STARRS 
and LSST discover more. This discussion has also focused on currently active MBCs. If these ob- 
jects are recently activated by collisions, there must also be a large population of icy asteroids that 
have not been recently activated. In the above analysis, all asteroids in a given target category are 
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assumed to be dormant MBCs, but in truth, the extent of this dormant population is unconstrained, 
except that it must assuredly be larger than the active population. In fact, if recent thermal models 
are correct, the population of dormant MBCs could, in principle, span the entire main belt (Fig. 4 
in ISchorghofer] [2008) 1. 



4.4. Impact Circumstances and Environment 



Hsieh et al. (2004 1 estimated a peak mass loss rate for 133P in 2002 on the order of ~0.01 kg s 
corresponding to an effective area of exposed water ice of up to a few hundred square meters 
(assuming a gas-to-dust ratio of 1), or an active "vent" (assumed here to be circular) with a radius 
of r,, ~ 10 m. Assuming that the vent was collisionally produced, this amount of exposed ice 
requires a crater of that size or larger. The size of the required crater is difficult to constrain as it 
depends on several factors, such as surface topology and the depth of the subsurface ice. If the layer 
of inert surface material is thick, volatile material may only become exposed in the deepest portion 
of the crater. Furthermore, a rapid expected mantling timescale (r„, ~ 15 days for a r ~ 2 km 



nucleus; Jewitt 1996 2002 1 suggests that much of the ice initially exposed by an activating impact 
is likely to be quickly reburied by rubble mantle growth, necessitating an initial crater that is larger 
than the calculated surface area of exposed ice needed to sustain the observed rate of sublimation. 

In order to calculate the size of the impactor required to create a reasonably-sized vent on an 
MBC, assuming that these vents are in fact collisionally-produced, one must first test whether the 
creation of 133P's active vent was likely to be gravity -dominated or strength-dominated using 



R 



3Y 



(4) 



where G is the gravitational constant, pa is the asteroid's mean bulk density (taken here to be pa - 
1300 kg m"^), is the asteroid radius, Y is the strength parameter (taken here to be F ~ 5 x 10^ Pa, 
a typical value used for porous silicates), and R > \ for gravity-dominated crater formation and 
< 1 for strength-dominated crater formation ( |Davis|[T999| l. The density used here is a minimum 
critical density required for 133P to be stable against centrifugal forces induced by its rapid rotation 



(Hsieh et al. 2004 1, but is also consistent with the density of pa - 1300 ± 300 kg m ^ found for 
(253) Mathilde (another C-type asteroid) from a flyby by the Near-Earth Asteroid Rendezvous 



(NEAR) spacecraft ( Veverka et al. 1999 1. One finds R ^ Ax 10"^ for the minimum inferred crater 
size of r,, ~ 10 m, and R x 7 x 10"^ for a crater as large as 133P itself {i.e., ~ 1.9 km). These 
results indicate that, due largely to 133P's extremely low surface gravity, any cratering on the object 
likely occurs in the strength-dominated regime. 

One therefore should use a crater scaling relation for the strength-dominated regime. 



C'd \ ^P' 



6pfl \ ' ( PiV, 



(5) 



where r, is the impactor radius, p,- ~ 2600 kg m ^ is the assumed impactor density, v,- ~ 5 km s ' 
is the assumed impactor velocity, and and j8 are material-dependent constants, adopted here to 



be = 1.5 and p — 0.165 (O'Brien et al. 2006 1. One thus finds that an impactor with a radius 
of at least r, ~ 0.6 m is needed to excavate an active site with a radius of r,, ~ 10 m. Using this 
radius as the minimum possible size for the active site, and taking this calculation as an order-of- 
magnitude estimate, hereafter I adopt r, ~ 1 m as the approximate size of the impactor responsible 
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for coUisionally activating 133P, and furthermore, assume this activating impactor size to also be 
typical for all other MBCs. 

Unfortunately, the population of meter-sized bodies in the asteroid belt is poorly constrained, 
meaning that impact rates by such small impactors are likewise poorly constrained. Extrapolations 
of the known main-belt population (which extends down to approximately kilometer-sized bodies), 
however, indicate that there may be approximately n, ~ 10'^ meter-sized bodies in the main belt, 
where both Cheng ( 20041 1 and Bottke et al. (2005i find an approximate sub-kilometer main-belt 
size distribution function of n,- ~ kr7^'^, where k = 3.2 x 10^. Approximating the main belt as a 
cylinder of r = 3.3 AU and h ~ 1.5 AU (where h is the approximate thickness of the asteroid belt 
measured perpendicular to the plane of the solar system), with an r = 2.1 AU cylindrical hole in 
its center, and assuming objects are distributed uniformly throughout the main belt, one derives an 
approximate impactor number density of A^,- ~ 10"^^ m"^. 

Then, using 

Ti ~ (v.^inrlNi)-' (6) 

for 133P with an effective radius of r„ - 1.9 km and having typical relative encounter velocities 
of Vrei ~ 5 km s ' with other main belt objects (e.g., Farinella & Davis 19921, one obtains an 
approximate impact time scale per object of t, ~ 2 x 10'^ s, or about one impact by a meter- 
sized projectile every 6000 yr (~1000 orbits). This should be regarded as an order-of-magnitude 
estimate only, as it is based on an assumption of uniform impactor density throughout the main 
belt (certainly untrue), which in turn is based on an extrapolation of the main-belt size-frequency 
distribution. For reference, following analogous lines of reasoning, order-of-magnitude timescales 
for impacts by 0.1 -meter-sized and 10-meter-sized projectiles are r, ~ 10 yr (~ 2 orbits) and 
Tj ~ 1 Myr (~ 2 X 10^ orbits), respectively. 



4.5. Active Lifetimes 

Like typical crater and impactor sizes for the MBCs, the typical lifetime of an active vent is difficult 
to constrain and is Ukely highly variable for individual situations. This is in fact the case for all 
comets (cf. [Sekaninaj 1 1 990| ), as numerous variables (e.g., location and depth of the activating crater, 
local ice content, local surface topography, orbital obliquity, perihelion distance, nucleus rotation 
rate) affect a particular vent's lifetime. Clearly a vent's active lifetime depends on the area of ice 
exposed, the size and morphology of the underlying volatile reservoir, and the depth of penetration 
of the activating impact. Surface topography near the active site and orbital parameters will affect 
the amount of solar radiation received, affecting the intensity and duration of sublimation activity. 
In the case of 133P, which has a rotation period of 3.471 hr ( jHsieh et a/.||2004| , centrifugal forces 
induced by rapid rotation could affect a vent's active lifetime by inhibiting rubble mantle growth, 
slowing this means of cometary deactivation. 

As directly calculating the lifetimes of active sites on MBCs is difficult to do with any confi- 
dence, an order-of-magnitude estimate can instead be derived from the statistics of the HTP survey 
itself. Earlier ( |4.3| l, it was estimated that 100 of the 8647 kilometer-scale, low-inclination, outer 
belt asteroids known at the time when the HTP was in progress could be currently active. Making 
the simplifying assumption that all asteroids in that population are equally likely to become active 
over their lifetimes, one finds that, on average, each individual asteroid is active over ~ 0.01 
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of its total lifetime. However, 133P is currently observed to be active over forh ~ 0.25 of its orbit. 
Thus, assuming 133P to be typical of all kilometer-scale, low-inclination, outer belt asteroids and 
if an activating collision occurs every ~1000 orbits, as estimated above, then the activity initiated 
by each collision should persist for ~40 orbits, or about 250 yr 

For MBCs with active orbit fractions of f orb - 0.1 or f„rb - 0.5, activity would instead be 
expected to persist for 600 yr and 120 yr, respectively. Additionally, if HTP survey statistics are 
incorrect by a factor of 2 (due, for example, to the extremely small-number statistics associated 
with a single detection, the uneven sensitivity of HTP observations, or other factors), the derived 
active vent lifetime could vary from 60 yr (using f,oi - 0.005 and forb - 0.5) to 1200 yr (using 
ftot - 0.02 and forb - 0.1). The typical lifetime of an active vent can therefore be estimated to be 
Tactv ~ 10^ - 10^ yr. 

If it is instead assumed that MBC activity is caused by 0.1 -meter-sized impactors, the higher 
impact rate of these more populous impactors implies Taav ~ 0.1 - 1 yr (following the line of 
reasoning above, using /,„, = 0.005 - 0.02 and /„,.;, = 0.1 - 0.5). If 10-meter-sized impactors are 
assumed to be the dominant activating population, the derived typical active lifetime is Tgav ~ 
10"^ - 10^ yr The active lifetime of 133P has already been observed to be at least 10 yr, well 
in excess of the lifetimes implied by 0.1 -meter-sized activating impactors. Meanwhile, the active 
lifetimes implied by 10-meter-sized activating impactors are implausibly comparable to or greater 



than the lifetimes estimated for classical comets from the outer solar system (~ 10 yr; e.g., Levison 



& Duncan 1997 1 that presumably contain much more ice, spend far smaller fractions of their orbit 



actually actively sublimating, and have only recently been perturbed into Sun-approaching orbits. 



Thus, the impactor size implied by the estimated crater size on 133P (§4.4 1 is consistent with 
impactor sizes implied by HTP statistics. 

4.6. Possible MBC Origin 

Having derived a self-consistent model of the activation and general properties of MBC activity, 
however, I turn to the issue of the origin of these bodies themselves. The surface area of 133P is 
approximately given by As - A-nr^ - 4.5 x 10^ m^. If a vent with r,, = 10 m is created every 
~6000 yr and the probability of impacts creating such vents is assumed to be equal over 133P's 
surface, one would expect the potentially volatile fraction (Fmi) of 133P's surface (i.e., the portion 
not yet devolatilized by impacts) to follow the exponential decay function 

F„„ = e-A' (7) 

where the fractional devolatilization rate, /a, can be approximated by /a - AsTT\nr-^.y\ and t 
is the elapsed time since 133P's formation. One then finds that 133P has a surface devolatization 
half-life of fi/2 ~ 700 Myr, and should have less than 10% of its surface unaffected by impacts 
after 2.5 Gyr, its current age if it is assumed to be an original fragment from the destruction of the 



Themis family parent body (§ 1.3 i 



If 133P is in fact a primordial member of the Themis family, the present-day timescale for 
impacts by meter-sized projectiles on still- volatile portions of its surface should be r, = Ti/Fmi ~ 
6 X lO'* yr. Assuming the activity detection statistics of the HTP are correct, this longer impact 
timescale then implies correspondingly longer active vent lifetimes of T„ctv ~ 10^-10^ yr. Such vent 
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lifetimes are not inconceivable but may be considered improbable given that they are approaching 



the estimated active lifetimes of much icier and better-preserved classical comets (c/. §4.5 i. 

Thus, for 133P to become activated on a plausible timescale without requiring implausibly long 
vent lifetimes, it may need to be a secondary (or even multi-generation) fragment from a recent 
disruption event. It and other MBCs could be sheared-off surface fragments or interior fragments 
of their respective parent bodies. In either case, the ice that they contain could have previously 
been located at great depths in those parent bodies, preserved against both solar heating and impact 
devolatilization, and now only recently found in shallow subsurface reservoirs. 

The conclusion that 133P could be a recently-produced fragment of a larger asteroid is sup- 
ported by [Nesvorny et af. (2008i who found that 133P may be a member of the newly-identified 
Beagle family (named for Themis-family asteroid (656) Beagle), which they believe formed 
< 50 Myr ago from the catastrophic disruption of a parent object with a diameter of D > 20 km. 
The other MBCs are not members of this new family nor are they known to be linked with any 
other known small families or clusters, but this does not necessarily preclude their membership in 



small families that have yet to be identified. Studies by [Cheng (2004 1 and Bottke et al. ( 2005 | l also 
predict that the coUisional lifetime against destruction of a D = 4 km asteroid should be ~2 Gyr, 
indicating that if 133? and 176P were formed in the initial Themis family fragmentation event, 
both should have been collisionally destroyed by the present day. 



5. Summary 

I describe the design, execution, and analysis of the Hawaii Trails Project survey that led to the 
identification of the class of main-belt comets. I note the following key points regarding the survey 
itself: 

- In an effort to find other 133P-like comets in the main asteroid belt, 657 observations were made 
of of 599 unique asteroids belonging to five major target categories which included the Themis, 
Koronis, and Veritas asteroid families, the Karin cluster (a subset of the Koronis family), and 
kilometer-scale, low-inclination, outer-belt asteroids, where the last category also includes all 
Themis family members and was formulated after the discovery of the second- and third-known 
MBCs to match the orbital characteristics of all three MBCs known at the time. 

- One surveyed object (176P; also asteroid 118401) was found to be active by this survey, and 
its discovery, in conjunction with the serendipitous discoveries of 133P and P/Read, led to the 
identification of the new class of main-belt comets. 

- There are strong observational selection effects, some intentional and some unintentional, that 
must be considered when interpreting the statistical imphcations of the HTP. Ultimately, the 
HTP should only be considered to be representative of the populations of objects actually sur- 
veyed. Even then, it should be noted that only one detection was made, and so the standard 
caveats associated with small number statistics apply. 

- The true abundance and distribution of main-belt comets in the main asteroid belt may not be 
truly known until the data from deep all-sky survey telescopes hke Pan-STARRS and LSST 
become available. These survey telescopes will likely find many more comets in the region 
of the asteroid belt where 133P, P/Read, and 176P are found, but may find comets elsewhere 
as well, such as where P/Garradd is found. Any MBCs discovered by these survey telescopes 
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will still require targeted foUow-up observations by human observers in order to confirm the 
presence of activity and determine physical and photometric properties. Dynamical analysis of 
newly discovered MBCs will also be essential for determining whether they are native to the 
regions where they are discovered or whether they may have migrated from elsewhere. 

I also make the following order-of -magnitude generalizations about the MBCs from the results 
of the HTP and currently available data on the currently best-characterized MBC, 133P (though the 
usual caveats associated with small-number statistics apply): 

- HTP survey statistics imply that there could be ~100 currently active MBCs among the 
kilometer-scale, low-inclination, outer belt asteroid population. Additionally, upper limits of 
150, 40, and 15 active objects among kilometer- scale members of the Koronis family, Karin 
cluster, and Veritas family, respectively, are found. The HTP provides no information about the 
prevalence of active MBCs outside these target categories. It also provides httle information 
about the size of the population of dormant MBCs {i.e., asteroids that contain subsurface ice 
but are currently inactive). 

- The estimated size of the active area on 133P is found to correspond to an impactor approxi- 
mately 1 m in size creating an active vent approximately 10 m in radius, assuming that the active 
area is collisionally created. This result is found to be consistent with the expected impact rate 
of projectiles of that size in the main belt and the cometary detection rate of the HTP. 

- Assuming that meter-sized impactors are responsible for initiating MBC activity, it is estimated 
that kilometer-scale, low-inchnation, outer-belt MBCs experience a potentially activating im- 
pact roughly every 6000 yr and remain active for approximately 10^ - lO' yr after each coUi- 
sional activation. 

- The estimated rate of impacts and estimated typical sizes of resulting active sites imply that 
133P should become significantly devolatilized over the course of 2.5 Gyr, suggesting that it 
may not be a primordial member of the ~ 2.5-Gyr-old Themis family. It is suggested instead 
that 133P, and possibly the other MBCs as well, are fragments of more recent disruptions. 
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Fig. 1. Distribution of observed objects by target category, where all Themis family members are 
also low-inclination outer-belt asteroids, and all Karin cluster members are also Koronis family 
members. Observations for which activity detection limits are fainter than 22.3 mag arcsec~\ 
25.5 mag arcsec~\ and 26.0 mag arcsec"^ are shown in Ught gray, medium gray, and dark gray, 
respectively. 
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Fig. 2. Semimajor axis versus eccentricity distribution of HTP survey targets, where proper orbital 
elements are plotted when available and osculating orbital elements are plotted for objects for 
which proper orbital elements are not available. Proper orbital elements for the background main- 
belt asteroid population are plotted using small black dots, while orbital elements for observed 
HTP targets are marked as labeled. Orbital elements of the currently-known MBCs are plotted for 
reference, although only 176P was observed as part of this survey. 
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Fig. 3. Semimajor axis versus inclination distribution of HTP survey targets, where proper orbital 
elements are plotted when available and osculating orbital elements are plotted for objects for 
which proper orbital elements are not available. Proper orbital elements for the background main- 
belt asteroid population are plotted using small black dots, while orbital elements for observed 
HTP targets are marked as labeled. Orbital elements of the currently-known MBCs are plotted for 
reference, although only 176P was observed as part of this survey. 
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Fig. 4. Orbital and physical properties of surveyed targets, including (a) semimajor axis (a) in 
AU, (b) eccentricity (e), (c) inclination (/) in degrees, and (d) absolute magnitude (//y). Proper 
orbital elements are used when available, while osculating orbital elements are used for objects 
for which proper elements are not available. Properties of the currently-known MBCs (where the 
absolute magnitude of P/Garradd is a lower limit) are marked for reference, although only 176P 
was observed as part of this survey. Observations for which activity detection limits are fainter than 
22.3 mag arcsec"', 25.5 mag arcsec"\ and 26.0 mag arcsec"^ are shown in light gray, medium 
gray, and dark gray, respectively. 
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Fig. 5. Observational circumstances of surveyed targets, including (a) heliocentric distance (R) 
in AU, (b) true anomaly (v) in degrees, (c) orbit plane angle (apf) in degrees, and (d) apparent 
/?-band magnitude. The observational circumstances of the currently-known MBCs during their 
observed active periods are marked for reference, although only 176P was observed as part of 
this survey. Observations for which activity detection limits are fainter than 22.3 mag arcsec"', 
25.5 mag arcsec ', and 26.0 mag arcsec"' are shown in Ught gray, medium gray, and dark gray, 
respectively. 
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Table 1. Telescope Summary 



Telescope 



Instrument 



Pix. Scalf " 



FO\ 



Nights A'oi 



CTIO 1.0 m 
CTIO 1.0 m 
Gemini 8.1 m 
Keck I 10 m 
Lulin 1.0 m 
SOAR 4.1 m 
Subaru 8.2 m 
UH 2.2 m 
UH 2.2 m 



Apogee512 

Y4KCam 

CMOS 

LRIS (Red) 

VA1300B 

SOI 

SuprimeCam 

Optic 

Tek2048 



a'47 
^'289 
a' 146 
a'210 
a'516 

a' 154 
a'2o 
a' 14 
a'2i9 



4:ox4;o 

19:6 X 19:6 
5:5x5:5 

7:2x7:2 

11:2 X 11:2 
5:3 X 5:3 
34' X 27' 
9:6x9:6 
7:5 X 7:5 



4 
1 
2 
7 
4 
1 
9 
52 



67 
24 
36 
52 
36 
28 
47 
48 
319 



" Pixel scale of instrument in arcsec pixel"' 
* Field of view in arcmin in the plane of the sky 
Total number of objects observed 
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Table 2. Observing Run Summary 



UTDate 



Telescope Instrument Weathei ' 



Moor 



Seeinj 



riob. 



2003 Sep 29 


UH 2.2m 


Tek2048 


clear 


N+3 


1.0 


3 


2003 Sep 30 


UH 2.2m 


Tek2048 


clear 


N-i-4 


1.1 


4 


2003 Dec 12 


UH 2.2m 


Tek2048 


clear 


N-11 


1.1 


3 


2003 Dec 13 


UH 2.2m 


Tek2048 


clear 


N-10 


1.5 


1 


2004 Apr 14 


UH 2.2m 


OPTIC 


clear 


N-5 


1.3 


2 


2004 Apr 15 


UH 2.2m 


OPTIC 


clear 


N-4 


1.5 


2 


2004 Jul 23 


UH 2.2m 


OPTIC 


clear 


N+6 


0.8 


1 


2004 Jul 24 


UH 2.2m 


OPTIC 


clear 


N-t-7 


0.9 


1 


2004 Jul 25 


UH 2.2m 


OPTIC 


clear 


N+8 


0.6 


4 


2004 Jul 26 


UH 2.2m 


OPTIC 


clear 


N-i-9 


0.8 


2 


2004 Aug 08 


UH 2.2m 


OPTIC 


clear 


N-8 


0.8 


13 


2004 Aug 09 


UH 2.2m 


OPTIC 


clear 


N-7 


0.6 


7 


2004 Aug 10 


UH 2.2m 


OPTIC 


clear 


N-6 


0.8 


16 


2004 Aug 13 


UH 2.2m 


Tek2048 


cirrus 


N-3 


1.2 


4 


2004 Oct 14 


UH 2.2m 


Tek2048 


cirrus 


N+0 


1.0 


8 


2004 Nov 04 


UH 2.2m 


Tek2048 


cirrus 


N-8 


1.3 


10 


2004 Nov 05 


UH 2.2m 


Tek2048 


cirrus 


N-7 


1.2 


6 


2004 Nov 06 


UH 2.2m 


Tek2048 


cirrus 


N-6 


1.1 


6 


2004 Nov 07 


UH 2.2m 


Tek2048 


cirrus 


N-5 


1.5 


3 


2005 Apr 10 


UH 2.2m 


Tek2048 


clear 


N+1 


0.7 


9 


2005 Apr 12 


UH 2.2m 


Tek2048 


clear 


N-H3 


0.9 


8 


2005 Apr 13 


UH 2.2m 


Tek2048 


cirrus 


N+4 


0.8 


17 


2005 Apr 17 


UH 2.2m 


Tek2048 


clear 


N+8 


1.1 


2 


2005 Apr 30-01 


CTIO 1.0m 


512x512 


clear 


N-8 


1.4 


4 


2005 May 01-02 


CTIO 1.0m 


512x512 


cirrus 


N-7 


1.2 


11 


2005 May 04-05 


CTIO 1.0m 


512x512 


cirrus 


N-4 


1.9 


7 


2005 May 06-07 


CTIO 1.0m 


512x512 


cirrus 


N-2 


1.7 


12 


2005 May 07-08 


CTIO 1.0m 


512x512 


cirrus 


N-1 


1.9 


1 


2005 May 08-09 


CTIO 1.0m 


512x512 


cirrus 


N+0 


1.6 


11 


2005 May 09-10 


CTIO 1.0m 


512x512 


cirrus 


N+1 


1.5 


9 


2005 May 10-11 


CTIO 1.0m 


512x512 


cirrus 


N+2 


2.0 


12 


2005 May 16 


UH 2.2m 


Tek2048 


clear 


N+8 


1.1 


5 


2005 May 27 


UH 2.2m 


Tek2048 


clear 


N-11 


0.9 


10 


2005 May 28 


UH 2.2m 


Tek2048 


clear 


N-10 


0.9 


10 


2005 May 29 


UH 2.2m 


Tek2048 


cirrus 


N-9 


0.7 


4 


2005 May 30 


UH 2.2m 


Tek2048 


cirrus 


N-8 


0.9 


11 


2005 May 31 


UH 2.2m 


Tek2048 


cirrus 


N-7 


0.9 


5 


2005 Jun 08 


Keck 10m 


LRIS 


cirrus 


N+1 


1.1 


19 


2005 Jul 06 


UH 2.2m 


Tek2048 


cirrus 


N+0 


1.1 


7 
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Table 2. continued. 



UT Date 


Telescope 


Instrument 


Weather 


Moon 


Seeing 


nobj 


2005 Jul 07 


UH 2.2m 


Tek2048 


clear 


N-Hl 


0.8 


6 


2005 Jul 09 


UH 2.2m 


Tek2048 


clear 


N+3 


1.7 


3 


2005 Jul 10 


UH 2.2m 


Tek2048 


clear 


N+4 


0.8 


7 


2005 Aug 25 


UH 2.2m 


Tek2048 


clear 


N-10 


0.8 


10 


2005 Aug 26 


UH2.2m 


Tek2048 


clear 


N-9 


0.6 


13 


2005 Aug 27 


UH2.2m 


Tek2048 


clear 


N-8 


0.7 


8 


2005 Aug 28 


UH2.2m 


Tek2048 


clear 


N-7 


0.9 


14 


2005 Aug 29 


UH2.2m 


Tek2048 


clear 


N-6 


0.7 


8 


2005 Aug 30 


UH2.2m 


Tek2048 


cirrus 


N-5 


0.9 


6 


2005 Sep 01 


UH2.2m 


Tek2048 


clear 


N-3 


0.9 


2 


2005 Sep 24 


UH 2.2m 


Tek2048 


clear 


N-9 


0.9 


10 


2005 Sep 26 


UH 2.2m 


Tek2048 


clear 


N-7 


1.2 


5 


2005 Sep 27-28 


CTIO 1.0m 


Y4KCam 


cirrus 


N-6 


1.7 


5 


2005 Sep 28-29 


CTIO 1.0m 


Y4KCam 


cirrus 


N-5 


2.5 


5 


2005 Sep 29-30 


CTIO 1.0m 


Y4KCam 


cirrus 


N-4 


2.3 


7 


2005 Sep 30-01 


CTIO 1.0m 


Y4KCam 


cirrus 


N-3 


1.5 


7 


2005 Oct 20 


Lulin 1.0m 


VA1300b 


cirrus 


N-13 


1.9 


5 


2005 Oct 22 


Lulin 1.0m 


VA1300b 


cirrus 


N-11 


1.7 


6 


2005 Oct 23 


Lulin 1.0m 


VA1300b 


cirrus 


N-10 


1.8 


5 


2005 Oct 24 


Lulin 1.0m 


VA1300b 


cirrus 


N-9 


1.5 


8 


2005 Oct 25 


Lulin 1.0m 


VA1300b 


cirrus 


N-8 


1.5 


6 


2005 Oct 26 


Lulin 1.0m 


VA1300b 


cirrus 


N-7 


1.4 


5 


2005 Oct 27 


Lulin 1.0m 


VA 1300b 


cirrus 


N-6 


1.5 


1 


2005 Nov 02 


UH 2.2m 


Tek2048 


cirrus 


N+0 


1.5 


2 


2005 Nov 03 


UH 2.2m 


Tek2048 


cirrus 


N+1 


1.5 


6 


2005 Nov 04 


UH 2.2m 


Tek2048 


cirrus 


N-H2 


1.0 


6 


2005 Nov 05 


UH 2.2m 


Tek2048 


cirrus 


N-H3 


0.6 


6 


2005 Nov 06 


UH2.2m 


Tek2048 


cirrus 


N-i-4 


0.8 


9 


2005 Nov 07 


UH2.2m 


Tek2048 


cirrus 


N-i-5 


0.7 


4 


2005 Nov 21 


UH2.2m 


Tek2048 


clear 


N-10 


0.7 


1 


2005 Nov 26 


Gemini 8m 


GMOS 


clear 


N-5 


0.6 


36 


2005 Dec 24 


UH2.2m 


Tek2048 


clear 


N-7 


1.0 


11 


2005 Dec 25 


UH2.2m 


Tek2048 


clear 


N-6 


0.9 


7 


2005 Dec 26 


UH 2.2m 


Tek2048 


clear 


N-5 


1.0 


8 


2005 Dec 27 


UH 2.2m 


Tek2048 


clear 


N^ 


0.8 


6 


2006 Feb 08 


UH 2.2m 


Tek2048 


clear 


Nh-10 


0.9 


2 


2006 Apr 22 


UH 2.2m 


Tek2048 


cirrus 


N-6 


1.1 


4 


2006 Apr 23 


UH 2.2m 


Tek2048 


cirrus 


N-5 


0.9 


9 


2006 Apr 28 


UH 2.2m 


Tek2048 


cirrus 


N-hO 


1.0 


1 


2006 Apr 30 


UH2.2m 


Tek2048 


cirrus 


N-i-2 


0.7 


2 
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Table 2. continued. 



UT Date 


Telescope 


Instrument 


Weather 


Moon 


Seeing 


' ^Ou J 


2006 May 22 


UH 2.2m 


Tek2048 


cirrus 


N-5 


0.8 


2 


2006 Jun 26 


Subaru 


SuprimeCam 


cirrus 


N+1 


0.6 


47 


2006 Jul 01 


UH 2.2m 


Tek2048 


cirrus 


N+6 


0.9 


3 


2006 Oct 17-18 


SOAR 


SOI 


cirrus 


N-4 


1.3 


2 


2006 Oct 18-19 


SOAR 


SOI 


clear 


N-3 


0.8 


7 


2006 Oct 19-20 


SOAR 


SOI 


clear 


N-2 


1.1 


4 


2006 Oct 20-21 


SOAR 


SOI 


clear 


N-1 


0.8 


15 


2006 Nov 11 


UH2.2m 


Tek2048 


cirrus 


N-10 


1.1 


2 


2007 Jan 27 


Keck 


LRIS 


clear 


N+8 


0.9 


33 



Total number of objects observed 657 

^ Weather conditions at the time of observations 
^ Lunar phase expressed in offset from new Moon ("N") in days 
^ Typical FWHM seeing in arcsec 
Number of objects observed 
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